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ABSTRACT 

We present the first high-resolution A''-body simulations of the fragmentation of 
dark matter filaments. Such fragmentation occurs in top-down scenarios of struc- 
ture formation, when the dark matter is warm instead of cold. In a previous paper 
(Knebe et al. 2002, hereafter Paper I), we showed that WDM differs from the standard 
Cold Dark Matter (CDM) mainly in the formation history and large-scale distribution 
of low-mass haloes, which form later and tend to be more clustered in WDM than in 
CDM universes, tracing more closely the filamentary structures of the cosmic web. 
Therefore, we focus our computational effort in this paper on one particular filament 
extracted from a WDM cosmological simulation and compare in detail its evolution 
to that of the same CDM filament. We find that the mass distribution of the halos 
forming via fragmentation within the filament is broadly peaked around a Jeans mass 
of a few 10^ Mq, corresponding to a gravitational instability of smooth regions with an 
overdensity contrast around 10 at these redshifts. Our results confirm that WDM fila- 
ments fragment and form gravitationally bound haloes in a top-down fashion, whereas 
CDM filaments are built bottom-up, thus demonstrating the impact of the nature of 
the dark matter on dwarf galaxy properties. 
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> 1 INTRODUCTION 
• I— I 



The Cold Dark Matter crisis on small scales is far from be- 
ing resolved. First, the highest resolution simulations avail- 
able still favour "cuspy" galaxy haloes (Power et al. 2003) 
with an inner slope for the density profile « —1.2, whereas 
high resolution observations of low surface brightness galax- 
ies are best fit by halo "cores" (de Block & Bosma 2002; 
Swaters et al. 2003). Second, the dearth of low luminos- 
ity galaxy satellites found A''-body simulations remains, al- 
though squelching of dwarfs due to the reionisation of the 
universe at high redshift may provide a solution (e.g. Ben- 
son et al. 2002; Somerville 2002). However, it is unclear that 
such a process is the definite solution, because the low-mass 
dark satellite haloes left over still interact gravitationally 
with the luminous galaxies, exchanging energy and angular 
momentum with them. Such interactions could have a sig- 
nificant impact on these galaxies' properties, heating their 
disks and reducing their sizes. Direct observational probing 
of dark matter halo substructure through anomalies in the 
flux ratio of multiple-image quasar lenses could in princi- 
ple determine whether the scenario previously suggested is 
the correct one (Metcalf & Madau 2001; Dalai & Kochanek 
2002), but there still are major difficulties to overcome in or- 



der to obtain reliable estimates of the amount of mass locked 
in dark substructures (Chen et al. 2003). 

An alternate mechanism to solve the satellite problem, 
which does not suffer from these side effects, is to reduce 
the power of dark matter fluctuations on small scales. This 
can be accomplished by changing the nature of the dark 
matter, (e.g. Kaplinghat, Knox, & Turner 2000; Spergel & 
Steinhardt 2000; Colin et al. 2000; Bode, Ostriker & Turok 
2001; Avila-Reese et al. 2001), by advocating an inflationary 
model with broken-scale invariance (Kamionkowski & Liddle 
2000) or by simply assuming an "unusual" dip in the primor- 
dial CDM power spectrum on small scales (Little, Knebe & 
Islam 2003). The latest constraints from polarisation mea- 
surements of the Cosmic Microwave Background Radiation 
by WMAP (Bennett et ah, 2003; Spergel et al., 2003; Kogut 
et al., 2003) seem to exclude such models, with the major 
caveat that the current theoretical understanding of the ob- 
jects that produce the first ionising photons is at best crude. 
As an example, a reionisation redshift of 11 (on the lowest 
side of the WMAP published 95 % confidence level) is still 
reached in a WDM model for which the warmon mass is 
about 1 keV, provided the objects responsible for photoion- 
ising the universe at this redshift emit on average between 
100 and 1000 times more ionising photons per baryon than 
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normal stars, as is possible with a top-heavy IMF dominated 
by accreting black holes (Schneider et al. 2002). Moreover, 
the constraints become less restrictive as the mass of the 
warmon increases. 

More sophisticated numerical simulations of liigh- 
redshift WDM-like universes, combined with a better un- 
derstanding of the reionising photon emitters are required 
in order to set definitive constraints on the nature of the 
dark matter from polarization measurements. This impor- 
tant point set aside, previous simulations (sec Paper I; Bode, 
Ostriker & Turok 2001; Avila-Reese et al. 2001) have indi- 
cated that deviation of WDM models from the CDM hier- 
archical structxirc formation scenario should be more pro- 
nounced in low density regions (e.g. filaments, voids). More 
specifically, in these simulations, most of the halos close to 
the resolution limit (i.e. less than 100 paxticles) seemed to 
be forming at smaller redshifts than in the standard CDM 
models, and in a top-down fashion through the break-up of 
filamentary structures. This phenomenon is of prime interest 
for galaxy formation and evolution models as it predicts that 
dwarfs properties strongly depend on the nature of the dark 
matter. It also gives rise to a series of interesting questions: 
how much do filaments fragment? When do they fragment 
the most? What is the mass distribution of the fragments? 
To tackle this issue we identified one of the filaments in the 
previous CDM and WDM low resolution runs presented in 
Paper I and resimulated it with a mass resolution 64 times 
higher than the original. 

This paper proceeds as follows: in Section 2 we briefly 
give the model parameters and describe the iV-body- 
simulations which we use in Section 3 to calculate the prop- 
erties of gravitationally bound halos found in the filament 
and of the filament itself. Finally, we present our conclusions 
in Section 4. 



2 AT-BODY SIMULATIONS 

The AWDM and the ACDM cosmological parameters used 
in this paper are the same as those of the AWDM2 and 
ACDM runs presented in Paper I, namely, fJo = 1/3, Ao = 
2/3, (78 = 0.88, h = 2/3, and mwoM = 0.5keV. We note 
that they are in good agreement with the latest results from 
WMAP (Bennett et al. 2003). 

As outlined in Paper I, initial conditions for our lower 
resolution simulations were generated by packing particles 
fi-om higher resolution initial conditions into groups. For the 
simulations presented in this paper, we simply identified a 
filament at 2: = in our 25 Mpc low resolution box, 
tracked back in time the Lagrangian volume its particles oc- 
cupied at 2 = 35, and unpacked the low-resolution particles 
into high-resolution particles in that region. In practice, our 
resimulation of the filament resulted in a 64-fold increase 
in mass resolution from nip ~ 7 x 10* M© (128^ parti- 
cles) in the low resolution runs to rup ~ 10^ Mq (512'^ 
particles) in the runs presented here. Each of our AWDM 
and ACDM filaments contains over 3 million high resolution 
particles and our force resolution is of the order i kpc 
in the densest regions. The total length of a filament can 
be represented by the dimensions xi, = 7, yi, = 4, zt, = G 
(in h^^ Mpc) of the rectangular box enclosing it. We took 
"snapshots" of the simulations at the following redshifts: 



Figure 1. Gray-scaled image of the logarithm of the density field 
around the filament at redshifts z = 1 (upper panel), 2 = 3 
(middle panel) and 2 = 5 (lower panel). 



z=5, 4.5, 4, 3.5, 3, 2.5, 2, 1.75, 1.5, 1.25, 1, 0.5, and 0. These 
re-simulations were undertaken using the publicly available 
Adaptive Mesh Refinement code MLAPM (Knebe, Green & 
Binney, 2001). 

Particle groups were identified using the Bound Density 
Maxima (BDM, Klypin & Holtzman 1997) method, in keep- 
ing with Paper I. The BDM code identifies local overdensity 
peaks by smoothing the density field on a particular scale of 
the order of the force resolution. We adopted a smoothing of 
10h~^ kpc for the runs under investigation. These pealcs axe 
prospective halo centres. For each of these halo centres we 
step out in (logarithmically spaced) radial bins until the den- 
sity reaches Phaio('"vir) = AvirPfc, where pi, is the background 
density and Avir ~ 340. This defines the outer radius rvir 
of the halo. Once we know the outer radius we are able to 
calculate internal properties of the halos. Note that during 
this group identification procedure we also check if each in- 
dividual particle truly belongs to the halo by comparing its 
velocity to the local value of the escape velocity. We exclude 
unbound particles from the haloes. This iterative method 
ensures that we only retain gravitationally bound objects. 
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Figure 2. Number of gravitationally bound halos located with 
the filament as a function of redshift. We employed a mass cut 
of M > lO^/i-i Mq (thick lines) and M > W^°h-^ Mg (thin 
lines), respectively. 

3 RESULTS 

3.1 Dark Matter Density Field 

We begin our investigation by looking at the large-scale 
structure density field within and around the filament. Fig. 1 
shows the particle distribution for both models with each 
particle being grey-scaled proportionally to the logarithm of 
the local (over-) density. A striking feature of Fig. 1 is the 
marked granularity of the density field in the CDM model 
at all redshifts in contrast with WDM. This efi^ect reflects 
the lack of filtering of power on small scales in CDM as 
opposed to WDM. However, while the WDM filament is ob- 
viously different from its CDM counterpart at redshift z=5, 
they are very similar at 2=1. This similarity is all the more 
remarkable that the logarithmic scale used to plot the den- 
sity fields shown in Fig. 1 artificially enhances the density 
contrast, thus highlighting differences between the two uni- 
verses that would have gone unnoticed had the scale been 
linear. 

A further difference of note is the appearance of individ- 
ual (low-mass) halos during the course of the WDM filament 
evolution, whereas dense objects only seem to be merging in 
the CDM filament, therefore decreasing steadily in number. 

We have not presented results at redshift z — because 
the picture at this epoch is qualitatively similar to that at 
reshift z = 1: most of the interesting processes leading to 
different looking WDM and CDM filaments occur at high 
redshifts {z > 1). 

3.2 Halo Abundance in the Filament 

To quantify the discussion in Section 3.1, we compute 
how the number of halos within the filament evolves 
with redshift. This number is plotted in Fig. 2, and con- 
tains all gravitationally bound particle groups above M > 
W^h~^ Mq (more than 100 particles - thick lines) and 
M > 10^"/i"^ Mg (more than 1000 particles thin lines), 
respectively. The former simply corresponds to the integral 
of the mass function over the entire halo mass range avail- 
able given our resolution. 
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Figure 3. Cumulative mass functions of all halos located within 
the filament (left panel) along with the mass function of 'new' 
halos as defined in the text (right panel). 



Fig. 2 demonstrates that independent of redshift, there 
are more halos in the CDM filament (as already expected 
from Fig. 1) than in the WDM filament. However, the ratio 
■^haVo^Z-^haio'^i whlch appears to be nearly independent of 
the mass cut employed, drops significantly from about 50 at 
z = 5 to approximately 7 at 2: = 2. This ratio is a factor two 
smaller for the higher mass-cut of 10^^ M©. The deple- 
tion of halos in both models from 2 = 2 to 2: = is readily 
explained when we consider that very few objects are cre- 
ated within the filament at low redshifts (see Fig. 5), and 
therefore, mergers must be the dominant process in deter- 
mining the number of halos in that redshift range. Indeed, 
as can be seen in Fig. 1 haloes stream along the web-like 
structures to eventually become subsumed within larger ob- 
jects (galaxy groups or clusters) located at the intersections 
of the filamentary network. 

Fig. 3 displays the cumulative mass function of all ha- 
los identified within the filament (left panel) along with the 
mass function of 'new' halos (right panel). The definition of 
a 'new' halo is given in Section 3.4 where Fig. 3 is analysed 
in more detail, but for the moment one can simply view 
them as haloes which did not have any progenitor at the 
previous output. At this point of the analysis, we stress that 
the cumulative mass function of all halos in the CDM fila- 
ment closely follows a power-law (with a slope close to -1) 
whereas in the WDM case one observes a deficit of low-mass 
objects, in agreement with the results presented in Paper I. 



© 0000 RAS, MNRAS 000, 000-000 



4 Knebe A., Devriendt J. E.G., Gibson B.K., Silk J. 




redshift z 

Figure 4. Fraction of mass within gravitationally bound haJos of 
mass M > 10^ Mq as a function of redshift. 



3.3 Mciss of the Filament 

We next addres the question of the fraction of the filament 
material which gets locked up in galactic halos and how 
much remains smoothly distributed in the filament. To an- 
swer this question, we calculate the ratio of particles be- 
longing to halos with M > 10^h~^ Mq (identified using the 
BDM code, cf. Fig. 2) to the total number of particles in 
the filament. This requires a definition for the boundary 
of the filament which follows from a previous definition for 
voids: Einasto et al. (1998) define void particles as having 
Pp < pb where pp is the local density at the particle posi- 
tion and Pb the cosmological background density. We adopt 
the same prescription and call those particles that follow 
Pp > Pb filament particles. The density at the particle po- 
sition is interpolated from the density field calculated on 
a regular 256^ grid; as we restrict ourselves to only high- 
resolution particles we are at all times confined to the fila- 
ment region (by definition). 

The results of this analysis are plotted in Fig. 4. This fig- 
ure highlights the overabundance of filament halos in CDM 
with respect to WDM, as the fraction of material in particle 
groups is at all times higher in the former than in the lat- 
ter. Only at z < 1 does the ratio Mhaios/Affliament become 
comparable in both models, which confirms the similarity of 
these structures previously noticed at z = 1 (cf. Fig. 1). We 
note however that from Fig. 1 alone, one might have claimed 
that most of the material should already be in haloes at 
z = 5 in the CDM filament, but Fig. 4 argues against it, 
as only about 10 % of the particles belong to haloes more 
massive than M > W^h~^ M©, for this model, at this red- 
shift. Of course, in the CDM case, this percentage is quite 
sensitive to the mass cut-off adopted, in the sense that if we 
had chosen to consider haloes composed of 50 particles or 
more, it would have increased to about 20 %. Nevertheless, 
this implies that the gravitational field within the filament 
is still dominated by the smoothly distributed dark matter 
component up to redshifts a ~ 1, both for the CDM and 
WDM filaments, and that accretion and merging are indeed 
important since the total number of halos decreases below 



z = 2 in both structures whereas the mass of each filament 
which is locked up within those haloes increases. 



3.4 New Halos in the Filament 

This leads us to the fundamental question of this study, 
namely how and when does the WDM filament fragment 
to form haloes? In other words, how important is this "top- 
down" fashion of forming haloes compared to the more tradi- 
tional hierarchical " bottom- up" scenarios of CDM-like mod- 
els? To tackle this issue, we constructed the merger tree of 
each individual halo between consecutive redshift outputs 
and tagged a halo as 'new' at a given redshift when no pro- 
genitor could be found at the previous output. The require- 
ment for a progenitor to be identified as such was that it 
had to contain a minimum of 100 particles or 1/10 of the 
actual halo particles depending on which of these two num- 
bers was the smallest. If such a progenitor weis found the 
halo was marked 'old'. Fig. 5 shows the fraction of 'new' ha- 
los formed per unit time as a function of redshift. As can be 
seen from this figure, the interesting redshift range for this 
analysis is from z = 4.5 to z = 2, as one observes a signifi- 
cant halo formation rate within the WDM filament during 
this time period (cf. Fig. 5). 

In order to obtain results not affected by spurious for- 
mation or evaporation of haloes between timesteps, we have 
set a lower mass limit of IQ^'h^^ M© (which translates into 
100 particles). This means that if a halo is less massive than 
this cut-off it has been neglected in our analysis. However, 
this does not imply that its progenitor(s) had to contain 
100 particles. The mass limit for the progenitor(s) is the 
minimum of 100 and 1/10 of the actual halo particles and 
hence can be as low as 10 particles (minimal mass for which 
a bound object can reasonably be detected) . We emphasize 
that in the extreme case where a 10 particle progenitor of 
a halo is identified, it only affects the 'tag' of the halo cur- 
rently under consideration, meaning that it is classified as 
'old' instead of 'new'. More specifically, if such a progenitor 
turns out to be unphysical (a mere collection of 10 parti- 
cles which happen to pass by one another quite slowly) it 
biases the results towards fewer new halos being created in 
the filament. In this sense, the curves plotted in Fig. 5 are 
lower limits on the true number of new haloes in the mass 
range \(fh~^ M© and above, as we have made sure that the 
halos marked as 'new' were not a part of any structure at 
the previous redshift. Fig. 5 alone strongly suggests that the 
hulk of the halo population is forming between z = 4.5 and 
z = 2.5 within the WDM filament whereas the same figure 
points at the majority of halos being already in place in the 
CDM structure by z = 4. To be more explicit, this figure 
shows that the WDM halo formation rate at z = 3 is still 
high enough that it could account for all halos present at 
this redshift, had it been constant in the past. But we also 
know from Fig. 2 that the maximum number of WDM halos 
is broadly peaked around z = 2, and from Fig. 5 that the 
halo formation rate is not constant but steeply increasing 
with redshift for z > 2. Furthermore, Fig. 5 and Fig. 7 tell 
us that the halo merger rate is negligible compared to the 
halo formation rate at z > 2. Therefore, we are forced to 
conclude that most halos were indeed 'born' in the WDM 
filament between z = 5 and z = 2. The same analysis, ap- 
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Table 1. Characteristics of the 749 particles that end up forming 
the gravitationally bound halo at 2: = 2 shown in Fig. 6. 
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Figure 5. Fraction of halos per Gyr appearing in the filaments 
and for which no progenitor could be found at the previous red- 
shift (new haloes). Note that the fraction becomes greater than 
unity at high redshift. This simply means that had the actual halo 
formation rate been the same at earlier times, more halos would 
have been formed than are currently present in the filament. 



plied to the CDM filament leads to the claim that most halos 
are in place before z — 5. 

This answers the question as to when the halos form 
out of the WDM filament, but Fig. 6 shows an example of 
how such a fragmentation occurs. The sequence of panels 
in that figure traces back in time all the particles found in 
a typical new halo of mass M « 7.5 X lO^h''^ Mq at red- 
shift z = 2 (i.e. a halo which was identified to be a 'new' 
gravitationally bound object at this redshift using the BDM 
algorithm). It is obvious from Fig. 6 that particles which end 
up within the virial radius of this halo at 2; = 2 are part of 
the filamentary structure until the very moment when the 
collapse takes place. The quantitative information obtained 
by tracing the particles pertaining to the final bound halo 
back in time and provided in Table 1 backs up this quali- 
tative analysis. In particular, we computed in Table 1 the 
averaged local density (in terms of the background density) 
at all 749 particle positions, as well as the mean extent of 
the region occupied by those particles. As a reference we 
also listed the virial overdensity value Avir for each of the 
redshifts under consideration in the last column of Table 1. 

To obtain a physical explanation for the fragmentation, 
we can calculate how the Jeans mass, A/j, varies with red- 
shift in the WDM filament prior to its fragmentation, as- 
suming its density field to be close to homogeneous on scales 
larger than the Jeans scale Aj. These masses and scales are 
given by the following formulae (e.g. Binney & Tremaine 
1987): 




1h-'Mpc 



Figure 6. Example of a halo fragmentation within the WDM 
filament. All 749 particles which end up in the virial radius of the 
bound halo at z = 2 are shown in black. The grey-shaded area 
represents the density field within a filament slice. The thickness 
of this slice is equal to the maximal distance between halo parti- 
cles at any time. 



Mj = 



A.J 



^5/2 



6G3/2 



Gpb{l + S) 



(1) 



(2) 



where p = (1 + S)pb is the density of the filament and ct„ 
its averaged velocity dispersion. In the previous equations, 
S stands for the (over-)density contrast within the filament 
with respect to the mean density of the background universe, 
pb, at a given redshift. For our cosmological model, Eq. (1) 
and Eq. (2) become 



Mj 



5.1 X 



10^(1 



-3/2 



(1 + 5) 



-1/2 



Xj = 0.13 X {i + zy''^\i + 5)-''^^ 



km/s 



km/s 
Mpc 



(3) 



(4) 



A rough estimate of the mass of the fragments can be 
calculated by assuming that the average density contrast of 
the filament, 5 « 10 is a slowly varying function of red- 
shift between z = 5 and z — 2, and that typical veloc- 
ity dispersions are of the order of 20 km/s*. We then get 
8 X IO^Mq < Mj < 2 X 10^"Mq , to be compared to the mass 
functions for new haloes given in Fig. 3 (right panel) at red- 
shifts 3.5, 2 and 1. We observe that the masses of the new 
WDM halo are actually close to these estimated values for 
the Jeans mass, thus supporting the assumption that they do 
indeed form through a fragmentation of the filament caused 

* This value is based on the averaged velocity dispersion of par- 
ticles measured in the simulated filaments at local overdensities 
of around 10. 
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by gravitational instability. Now, the same range of values 
for 5, z and cr„ yields 53 kpc < Aj < 150 kpc which validates 
our assumption for the homogeneity of the WDM filament 
density field on these scales since the free-streaming length 
7?/ of a 0.5 kcV particle is i?/ ~ 300 kpc corresponding to a 
free streaming mass Mj ~ 4 x lO^'^M©. Note that this is not 
the case for the CDM filament as the free-streaming scale 
associated with the cold particle is much smaller than the 
Jeans length, and therefore we cannot consider the density 
field as homogeneous on these scales. Moreover, one consid- 
eration regarding the simulation technique needs to be borne 
in mind: our Af-body code MLAPM achieves high spatial reso- 
lution primarily in high density regions. This means that the 
force resolution can smooth out the density field on scales 
above the Jeans mass in lower density regions, with the re- 
sult that one would then tend to identify haloes as 'new' that 
should have already collapsed on smaller scales at earlier 
times. Nevertheless this possible numerical fragmentation is 
not too much of a concern in the present study, because 
even if one was to consider all the new haloes identified in 
the CDM filament at all redshifts as numerical artifacts, one 
still would have to invoke physical fragmentation to explain 
how the more massive new WDM haloes formed (see Fig. 3) . 



3.5 Merger History in Filament 

Parallel to fragmentation and creation of new haloes, we also 
measure the rate of major mergers within the filaments in 
both models. A major merger event is defined by the ratio 
of the two most massive progenitors: if this ratio is smaller 
than 3:1, then the merger is major, otherwise it is minor. As 
before, a lower mass cut-off of VfhT^ M© is applied for each 
halo. Fig. 7 indicates that at early times (high redshift) there 
is no major merger activity in the WDM filament whereas 
at later times the relative fraction of major mergers exceeds 
that found in the CDM scenario. This supports our previ- 
ous claim that haloes at high redshift [z > 3) do not form 
according to the traditional bottom-up hierarchy (i.e. via 
mergers) in the WDM filament. Moreover, since the mass 
distribution of the fragments formed through the top-down 
mechanism between z = 5 and z = 2 is broadly peaked 
around the Jeans mass (see Fig. 3), one naturally expects 
the major merger activity to reach a maximum around z = 2 
in WDM because this epoch combines a high number density 
of low mass halos with a small mass spread in their distribu- 
tion. In other words, at redshift higher than 2, quite a large 
fraction of small halos are still forming from fragmentation, 
and at lower redshift most of the formation and merging 
have already taken place, broadening the mass distribution 
function. As a result of these cfi'ects, the major merger frac- 
tion is strongly suppressed both at low and high redshifts. 
On the other hand, in the CDM filament, the fraction of 
major mergers is almost constant in time and reaches more 
modest values. This is the result of most of the haloes being 
already in place at high redshift {z > 5) which guarantees 
that the mass distribution of haloes will broaden smoothly 
over time as merging and accretion proceed. 



0.20 F 




no I \ 

12 3 4 

redshift z 



Figure 7. Fraction of merger events per Gyr for which the ratio 
of the two most massive progenitors is higher than 3:1 (major 
mergers) . 



4 DISCUSSION AND CONCLUSIONS 

We presented the first high resolution iV-body simulation of 

the top-down formation of dark matter halos through the 
fragmentation of a dark matter filament. This loads to a 
number of interesting results: the epoch of formation of small 
haloes is pushed back in time to lower redshifts (2 < 2; < 5) 
with respect to the standard CDM model, and the mass 
distribution of the fragments is broadly peaked around a 
Joans mass of a few 10^ Mq, corresponding to a gravita- 
tional instability of smooth regions with an overdensity con- 
trast around 10 at these redshifts. We need to undertake 
detailed hydrodynamic and star formation simulations to 
properly assess the impact of such a mode of formation on 
dwarf galaxies, but we anticipate the predictions in terms 
of clustering and the mean age of the stellar populations to 
be quite different from those of the standard CDM model. 
We defer this study to a companion paper. Reionisation of 
filamentary structures might not be a problem since we have 
shown that massive objects are already in place at z = 5. Of 
course, the higher the redshift the more difficult it will be 
to reionise the universe unless the decrease in the number 
of objects is compensated by an increase in the efficiency 
of the ionising sources. Our numerical simulations strongly 
suggest that spurious numerical fragmentation may occur in 
CDM models due to the softening of the gravity force. How- 
ever, this is not to be confused with the real fragmentation 
observed in the WDM simulation presented here because in 
AMR codes such as MLAPM the force resolution is a func- 
tion of the mass resolution. Therefore, two haloes with the 
same number of particles in the ACDM and AWDM sim- 
ulations are resolved with the same force resolution. This 
implies that if we were dominated by artificial fragmenta- 
tion in both ACDM and AWDM we should not observe 
different upper mass cut-offs in the mass distribution func- 
tions of new halos. Looking at Fig. 3 we clearly see that 
this is not the case: the mass distribution of new halos in 
ACDM drops off sharply at « 2 x 10^ M© whereas it 
extends to masses > W^°h~^ Mq in AWDM. Furthermore, 
the resolution of our new WDM haloes is high enough that 
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they should have been detected by the BDM algorithm at 
higher redshift had they been present at that time, since 
CDM haloes as massive are indeed classified as old ! Fi- 
nally, we point out that Melott (1982) already showed that 
galactic haloes wore forming via real fragmentation in Hot 
Dark Matter simulations where the resolution was several 
orders of magnitudes below what we achieve in this study. 
In light of these remarks, we conclude that while our results 
for the formation of massive (at least several hundred parti- 
cles) WDM haloes through fragmentation arc robust, a more 
detailed study is needed in order to really quantify numeri- 
cal fragmentation effects on smaller scales, as they might be 
relevant for studies of dwarf galaxies in cosmological con- 
texts where small haloes arc not sufficiently resolved. More 
specifically, wo believe these numerical effects will lead to a 
later formation of low mass haloes in moderately overdense 
regions. 
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